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Abstract. This paper is concerned with mathematical modeling of intelligent 
systems, such as human crowds and animal groups. In particular, the focus 
is on the emergence of different self-organized patterns from non-locality and 
anisotropy of the interactions among individuals. A mathematical technique 
by time-evolving measures is introduced to deal with both macroscopic and mi- 
croscopic scales within a unified modeling framework. Then self-organization 
issues are investigated and numerically reproduced at the proper scale, accord- 
ing to the kind of agents under consideration. 



1. Self-organization in many- particle systems 

One of the most outstanding expressions of intelligence in nonclassical physical 
systems, such as human crowds or animal groups, is their self-organization ability. 
Self-organization means that the individuals composing the system can give rise to 
complex patterns without using intercommunication as an essential mechanism. 

For instance, in normal conditions pedestrians are known to arrange in specific 
patterns, chiefly lanes (cf. Fig. [iji-b), as demonstrated by many experimental 
investigations [THl [201 1231 12H1 121] • Lane formation may be fostered by a suitable 
setup of the space, as reported in (THl 123]: a test performed in a tunnel connecting 
two subway stations in Budapest showed that a series of columns, placed in the 
middle of the walkway, induce pedestrians to organize in two oppositely walking 
lanes, preventing each of them to expand up to the full width of the corridor. More 
in general, lanes form also spontaneously, i.e., without the need for being triggered 
by environmental factors, provided the density of pedestrians is sufficiently large 
[2U] . This is particularly evident if one considers the case of two groups of people, 
walking in opposite directions, which meet and cross (see also [3S]). 

Grouping and self-organization are well known and largely observed also in ani- 
mals, see for example |34j . These phenomena are in fact ubiquitous, ranging from 
bird flocks in the sky to migrating lobsters on the sea floor. Many papers on this 
subject (see, among others, [71[II1[T1[I1[IS1[IZ1E11[3S1[M1I1H1EI|) proved by 
means of numerical simulations that few simple rules adopted by each animal can 
give rise to a complex organization of the whole group. Patterns commonly seen 
in nature are (see [26] for the nomenclature): globular clusters (e.g., starlings. Fig. 
[ij;, surf scoters. Fig. [iji), extended and front clusters (e.g., wildebeests, antelopes, 
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(b) (e) 

Figure 1. Self-organizations, (a), (b): oppositely walking lanes 
of pedestrians, (c): three-dimensional globular cluster of starlings. 

(d) : two-dimensional crystal-like globular cluster of surf scoters. 

(e) : a line of migrating African elephants. Reproduction of these 
pictures with kind permission of the respective copyright holders. 
Credits are in the acknowledgments at the end of the paper. 



pigeons), lines (e.g., penguins, lobsters, elephants, Fig. [T^), Vees, Jays, and ech- 
elons (e.g., geese). Motivations and convenience of these patterns are still under 
investigations, especially for Vees. 

It is commonly agreed that self-organization is the result of elementary actions 
that each subject performs to fulfill specific wills. Concerning pedestrians, the 
following basic guidelines can be identified: 

• The will to reach specific targets, e.g., an exit or a meeting point, which 
drives pedestrians along preferential paths, determined mainly by the ge- 
ometry and the spatial arrangement of the walking area. Unlike animals, 
pedestrians experience strong interactions with the environment, because 
they usually move in highly structured spaces scattered with all sorts of 
obstacles. 
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• The will to not stay too close to one another, with a preference for un- 
crowded areas (repulsion from other individuals). Pedestrians may agree 
to deviate from their preferred path, looking for free surrounding room. 

In addition, it is reasonable to believe that occasionally also a mild form of cohesion 
occurs, which translates the tendency of pedestrians to not remain isolated. This 
happens, for instance, in those groups whose individuals share specific relationships, 
such as groups of tourists in guided tours. 

For animals, the basic guidelines can instead be outlined as follows: 

• The will to not stay too close to one another, in order to avoid collisions 
{repulsion from other individuals). 

• The will to not remain isolated (group cohesion). Grouping is in fact ad- 
vantageous for many reasons, such as predator avoidance or food search, 
see [31]. 

Like pedestrians, also animals may want to reach some specific destinations (e.g., 
when they migrate). However, the direction of motion toward targets is usually 
almost constant for quite a long time, because animals move in large and basically 
obstacle-free environments. Therefore, this aspect plays a minor role in the descrip- 
tion of their self-organization, corresponding to a simple translation of the center 
of mass of the whole group. 

Self-organization can be broken under particular circumstances entailing a dra- 
matic change in the basic interaction rules discussed above. For pedestrians, an 
illuminating example is panic, when individuals tend to cram toward a common 
target (e.g., an escape) instead of seeking the less congested paths. Lane forma- 
tion is then ruled out and less organized patterns emerge in the crowd, probably 
in consequence of a strong simplification of the interaction rules. Nevertheless, the 
arrangement of the environment may help restore, at least partially, the normal 
order of the fiow. A typical example is that of an obstacle placed in front of an 
exit, which in some situations (presumably under panic conditions) may improve 
the flow of people through the exit itself, provided shape, dimensions, and position 
of the obstacle are accurately studied. This is a variant of the so-called Braess' 
paradox [31] , which states that a condition intuitively expected to lead to a worse 
situation may instead give rise to better outcomes. In animals, the changes in the 
behavior are even more evident. External conditions (e.g., presence of predators, 
weather), group tasks (feeding, exploring) or group speed can modify the interac- 
tion rules, leading to great modifications in the resulting patterns of the group. 
The environment usually does not interfere much, and pattern formation is mainly 
due to interactions among the individuals. 

In this paper we are concerned with mathematical modeling of the above-discussed 
systems, with the specific aim of describing the spontaneous emergence of self- 
organization. In particular, we focus on two basic characteristics of the interactions 
among the individuals, namely non-locality in space and anisotropy. The inclusion 
of these factors makes our models able to explain the differences observed in self- 
organization and pattern formation of various groups of agents in terms of different 
visual fields and sensing zones. 

At the same time, we introduce a modeling framework based on the measure the- 
ory, that allows for a unified formulation of macroscopic and microscopic models in 
any space dimension, and for a convenient numerics. This enables us to investigate 
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both macroscopic self-organization, typical of large crowds of pedestrians, and mi- 
croscopic self-organization, more specific of animal groups, using common modeling 
principles and tools. 

In more detail, the paper is organized as follows. After this introduction. Sect. [2] 
proposes a comparison between classical and intelligent particles, and highlights the 
main differences of the latter with respect to more standard systems dealt with by 
classical physics. In view of these considerations. Sect. [Sjdevelops some preliminary 
modeling strategies to enhance intelligence as a distinctive feature of the systems 
at hand. These serve as guidelines for Sect. [4] where the modeling technique 
by time-evolving measures is introduced and specific macroscopic and microscopic 
models are detailed. Numerical tests on the ability of these models to reproduce 
spontaneous self-organizing patterns for both human crowds and animal groups are 
performed and commented in Sect. [5j Finally, Sect. [6] draws some conclusions and 
briefly sketches research perspectives. 

2. Classical vs. intelligent particles 

A plethora of physical systems can be basically described as systems of inter- 
acting particles, think for instance of fluids, gases, and similar matters. Actually, 
also human crowds and animal groups are susceptible to this rough characteriza- 
tion, although it should be clear from the previous discussion that the particle 
analogy is only formal. Indeed, intelligence is what really makes the difference in 
this context: it gives rise to a decision-based dynamics determined by individual 
behavioral rules, rather than to a classical dynamics passively ruled by inertia. For 
this reason, human beings and animals are more properly characterized as Intelli- 
gent Particles (denoted IPs in the sequel for brevity). IPs can act directly on the 
system, rather than being passively subjected to the evolution itself, whence their 
ability to self-organize and to generate complex patterns. 

In the following, we summarize the main differences between classical and intel- 
ligent particles, in order to outline the main novelties posed by intelligent systems 
with respect to more standard frameworks. 

Robustness vs. Fragility. Classical particles are robust, in the sense that they in- 
teract almost exclusively through collisions. For instance, gas particles change 
direction of motion and velocity only when hit by other particles, or possibly when 
they collide with the walls of the container in which the gas is stored. Conversely, 
IPs are fragile, they try as much as possible to avoid mutual collisions as well as to 
steer clear of walls and obstacles scattered along their path. 

Blindness & Inertia vs. Vision & Decision. Classical particles are blind, indeed they 
have no information on the environment and on the distribution of the surrounding 
particles. Therefore, their dynamics is essentially ruled by inertia, i.e., a passive 
response to the mechanical cues coming from the exterior. On the contrary, human 
beings and animals feature specific visual fields, hence they can obtain information 
on the surrounding environment (e.g., presence of obstacles, of walls) and on the 
current distribution of (a certain number of) other agents. This information is then 
used to make decisions on the future individual evolution, generally by following 
some ordinary behavioral rules proper of the kind of agent under consideration. 
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Local vs. Nonlocal Interactions. Collisions among, say, gas particles require the 
latter to be sufficiently close to hit each other. Due to the typical size of gas 
molecules, by far much smaller than the environments where they flow, this has been 
classically understood as if the colliding particles were occupying the same spatial 
position at the moment of the impact (see e.g., [49 and the references therein), 
hence assuming the conceptual approximation of local, i.e., pointwise, interactions. 
Conversely, IPs do not interact mechanically, rather they are influenced by the 
presence of other individuals or objects a certain distance away, that they want 
either to approach or to avoid. The resulting interaction is thus nonlocal, because 
the agents do not need to be in contact to interact. 

We notice that there exist basically two types of nonlocal interaction: a metric 
one, such that each agent is influenced by all other agents located at a distance less 
than a given threshold; and a topological one, such that each agent is influenced by a 
given number of other agents, no matter how far they ar^ Topological interactions 
have been used in several models, e.g., HH [5T] , this idea being supported 

by experimental investigations like, for example, [2] on fish schools. Other models 
[TTJ [m [171 ES] do not include this feature in favor of a purely metric approach. 
However, recent results [3J show that starlings interact topologically with other 
six/seven group mates. It is reasonable to believe that topological interactions are 
common in nature, even in pedestrians, since they are mainly due to a limited 
capacity in processing the information. On the other hand, a purely topological 
interaction in a wide domain is not realistic, for an individual may not be concerned 
with very far mates (e.g., because it does not see them at all). 

Isotropy vs. Anisotropy. Another striking difference concerning the way in which 
classical and intelligent particles interact is that the first are isotropic, meaning that 
they are equally affected by mechanical cues coming from all directions. IPs are 
instead anisotropic, i.e., they are sensitive to stimuli coming from specific directions. 
This is partly related to the width of their visual field, which in pedestrians coincides 
with the half-space in front of them, whereas in animals covers often almost all the 
surrounding space. More precisely, the visual field must be intended as an upper 
bound for the regions in which cohesion and repulsion are active. Indeed, repulsion 
can be expected to be mainly felt against the IPs ahead rather than against those on 
the side or behind, especially in running people or fast-moving animals. Cohesion 
instead is active toward the group mates in front, if the goal is to follow the head 
of the group, or in every direction, if the goal is the unity of the group. 

The ideas of limited visual field, expressed in terms of blind rear zone, and of 
anisotropic sensing are quite common in the biological literature, see for instance 
[TU [TTJ H5j 1321 ESI [3S1 137] • However, they mainly play the role of passive features 
of the agents, and are not regarded as active features able to influence the shape of 
the group. 

Energy & Entropy vs. Self-organization. Particle collisions allow to describe the 
mechanics of classical systems by invoking the balance of linear momentum, as well 
as energy and entropy principles. A straightforward extension of these ideas to 
intelligent systems is not possible, because the ability of IPs to make subjective 



It is worth noting that in the biological literature authors often regard both metric and 
topological interactions as local models, in contrast to nonlocal models in which every agent 
interacts with all other group mates. 
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decisions continuously puts and removes energy in and from the system, in hardly 
quantifiable amounts. Entropy criteria may in turn be questioned, because entropy 
is classically related to the equiprobability of the states of a system, whereas self- 
organization promotes specific configurations of the IPs. The idea that, in living 
matter, energy and entropy play a somehow atypical role was already suggested by 
Schroedinger |46j in the Sixties, who formalized it through the concept of negative 
entropy. 

It is useful, at this point, to consider a very borderline system with respect to 
the "classical vs. intelligent" dichotomy, namely vehicular traffic. Traffic is, in 
principle, an intelligent mechanical system, because the intrinsic mechanics of the 
vehicles is tempered by the presence of the drivers, who determine a non strictly 
mechanical behavior of the cars. However, many successful models, relying on the 
fluid dynamical analogy of the flow of vehicles along a road, have been proposed in 
the literature (see \27\ I45j for comprehensive reviews on microscopic, kinetic, and 
macroscopic models). This approach has been possible, despite the nonstandard 
nature of the system at hand, because traffic is essentially one-dimensional, so that 
vehicles regulate their speed only. Such a dimensional constraint reduces signif- 
icantly the possibility of self-organization and pattern formation, therefore fluid 
dynamical modeling, entropy reasonings, and finally the development of models 
and theories based on nonlinear conservation laws are feasible. 

3. What mechanics for intelligent systems? 

The discussion of the previous sections has highlighted a crucial feature of in- 
telligent systems: their dynamics is only partially determined by classical (pseudo- 
)mechanical cues, therefore, in principle, it is not fully describable within the clas- 
sical Newtonian framework of point mechanics. Concerning this, we note that in 
the systems we are considering impulsive forces are often present. For example, an 
animal which starts moving reaches in a very short time its final velocity, which 
then remains constant for a while. Including in the models such impulsive forces 
through a Newtonian approach is quite difficult and, after all, needless. 

The dynamics of intelligent systems calls for a paradigm which enhances intel- 
ligence as a distinctive feature, rather than recovering it as a by-product of other 
principles. In this respect, let us consider that, from the mechanical/dynamical 
point of view, intelligence might be regarded as the ability of the agents to control 
their velocity actively, i.e., without the need for inertial accelerations caused by 
external actions. In order to take this aspect into account, we suggest to split the 
velocity v of the agents in two basic contributions: 

(1) V = w + ly, 

to each of which there correspond specific modeling strategies. In more detail: 

• w is the external velocity, i.e., the component of the total velocity depending 
essentially on actions exerted on the agents by the external environment. 
For example, in human crowds these may be the repulsion produced by 
obstacles and walls or the attraction produced by targets (e.g., doors, dis- 
plays, other people). In animal groups this effect is instead greatly reduced, 
because animals move in much less structured and generally clear environ- 
ments. For the external velocity, a Newtonian-like modeling is feasible, 
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hence w can be deduced from inertial reasonings up to a careful identifica- 
tion of the external actions. 
• v is the intelligent velocity, i.e., the component of the total velocity that 
the individuals control actively. The intelligent velocity is determined by 
the behavioral rules that the agents comply to, and by the decisions that 
they consequently make, therefore it need not be Newtonian. For instance, 
pedestrians and animals set up the intelligent velocity according to the 
occupancy of the surrounding space, so as to steer clear of congested areas 
while possibly preserving the compactness of the group. 
A convenient manner to model v is by an equation of state of the form 

Ht, x)^m]it, x), 

where t, x are time and space variables, respectively, Q comprises all state variables 
which can contribute to the determination of v (e.g., the distribution of the agents 
in a suitable neighborhood, the metric and topological structure of the sensing 
zone), and / is a functional relationship. 

Many microscopic models of animal groups available in literature actually fit 
into such framework, and the same is true for macroscopic models using hyperbolic 
partial differential equations, see e.g., [15l |47j and references therein. The forces 
do not appear exphcitly, the velocity being expressed directly as a function of the 
density of the animals, usually in a nonlocal way. On the other hand, different 
models at both the microscopic and the kinetic scale rely instead, at least formally, 
on a more Newtonian-like framework, in the sense that they recover the dynamics 
of the system from generalized forces responsible for the acceleration of the agents 
[51 [THl m] . Similar considerations apply also to models of human crowds, for an 
overview of which we refer to [^T]. In particular, macroscopic models derived by 
closing the mass conservation equation with suitable relations for the velocity fol- 
low the ideas outhned above, see e.g., El UHl [Ml ED- Other models resort 
instead more heavily to classical fluid dynamical analogies, indeed they supplement 
the mass conservation equation by a momentum balance equation invoking con- 
cepts like acceleration or generalized forces and pressures In these models, it 
is in principle possible to introduce topological interactions among IPs, as well as 
anisotropic sensing zones. Unfortunately, models based on PDEs and nonlinear 
conservation laws suffer from an important drawback: their mathematical treat- 
ment and numerical implementation get rather complicated in dimensions greater 
than one. Indeed, in realistic cases of interest for applications one must deal with 
possibly punctured two-dimensional domains (the holes representing the obstacles) , 
which requires to handle boundary conditions in a way that may not be immediately 
suited for hyperbolic equations. Therefore, we prefer to use an alternative model- 
ing strategy, which, at the same time, enables us to develop a modeling framework 
capable to treat both the macroscopic and the microscopic scale within a unified 
theoretical formulation. Then we will show the ability of specific models, deduced 
from this framework, to describe the emergence of self-organizing behaviors at the 
proper scale, according to the kind of agents under consideration. 

4. Mathematical modeling by time-evolving measures 

In this section we set up a modeling framework for the study of the behavior of 
IPs. In doing this, we resort to a discrete-time dynamics and describe the occupancy 
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of the space by the agents at different times via a sequence of positive Radon 
measures. This approach is inspired by [5], where the authors address rendez-vous 
problems for multi-agent systems, and has already been investigated in recent works 
|43l I44j . Here we show furthermore that also microscopic models for animal groups 
can be recast in the time-evolving measures framework. Therefore, we ultimately 
provide a unified modeling procedure to treat both macroscopic and microscopic 
intelligent systems. 

Let us consider a domain 51 C M"^, which represents the area where IPs are 
located and move. In our discussion, d is the dimension of the domain, from the 
physical point of view it may be c? = 1, 2, 3. A great advantage of our approach 
is that there are basically no differences in the theory for different dimensions. At 
every time instant n > 0, we define a Radon positive measure /i„ over the Borel a- 
algebra B{^), such that, for each measurable set E E B{il), the number fJ,n{E) > 
measures the occupancy of the area E by the IPs. In other words, the mapping 
fi„ provides the localization, i.e., the distribution, of the agents in the domain 
at time n. The evolution to the next time n + I depends on the dynamics of the 
system, that we describe via a motion mapping 7„ : il — s- il: 

(2) Jn{x) — X + Vn{x)At, X E ft. 

In practice, in the time step n ^ n + I the point x is advected to 7n(a^) by the 
velocity field u„ : ^ R'^. The duration of the time step is At > 0. Then, the 
measure /i„+i is constructed by pushing /i„ forward with 7„, that is /i„+i — jni^lJ'n 
or, more explicitly, 

(3) fin+i{E)^finh-\E)), VEeBin). 

This corresponds to the simple idea that the number of IPs contained in E at time 
n + 1 coincides with their number at the previous time n in the pre-image j~^(E), 
therefore Eq. (|3| expresses the conservation of the mass of IPs. By rewriting it in 
the equivalent form 

- M„(i?) = -[m„(7-1(£;^) nE)~ ^ir.{E- n i-\e))], 

we recognize the formal statement of a conservation law: the left-hand side repre- 
sents the variation of the measure of ii' in a single time step, and the right-hand side 
the difference between the outgoing flux, i.e., the measure of the set of points be- 
longing to E and mapped outside _B by 7„, and the incoming flux, i.e., the measure 
of the set of points not belonging to E and mapped inside £' by 7„. 

Mathematical models are deduced from the structure ([2])-([3| by prescribing the 
initial distribution of the IPs, that is the measure /ip, and by specifying the form 
of the velocity field w„. According to the discussion of Sect. [3] we represent as 
the superposition of external and intelligent contributions: 

Vn{x) ^W{x) + Vn[t^n]{x). 

The functional dependence of the intelligent velocity on the distribution of the 
agents, formally expressed by the square brackets in this formula, enables us to 
account for all of the factors having to do with the perception of the occupancy 
of the surrounding space from IPs (cf. Sect. [1. In more detail, we distinguish 
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Figure 2. Definition of the zone of cohesion by means of the ve- 
locity w, the topological radius Rc, and the span angle etc G (0, 27r]. 
Analogous considerations hold also for the zone of repulsion, up to 
considering a metric radius Rr and an angle a^. 



cohesive and repulsive effects in z^„, that we model as 

J \y — x\ 

for coefficients Fc > 0, Fr < 0. In this equation, Bc{x), Br{x) C are the zone 
of cohesion and the zone of repulsion of the agent x, respectively, which will be 
defined in the following. 

The strength of the cohesive velocity v'^ is chosen to be growing linearly with 
the distance among IPs within the zone of cohesion Bc{x). Cohesion is mainly 
topological, i.e., it involves a predefined number of IPs. By consequence, for a 
fixed positive p, the zone of cohesion is chosen to satisfy ^niBc{x)) < p, but its size 
may vary from point to point according to the distribution of the agents. However, 
a maximal size of Bc{x) exists, say s, beyond which interactions among IPs are 
inhibited due to the distance (cf. Sect. [2]), so that also C'^{Bc{x)) < s {C^ being 
the Lebesgue measure on M'*). In practice, Bc{x) is adjusted dynamically under 
the constraints 

(5) iin{Bc{x))<p and C\B,{x)) < s. 

Let us choose, for convenience, Bc{x) to be a circular sector of the ball of radius 
Rc > (to be detemined), centered in x and with central angle ac S (0, 27r] (Fig. 
[2]). Other choices are possible, for instance in [35] the authors use elliptical regions 
to take the body form of the agents into account. Due to the symmetry of the 
body, we assume that the sensing domain of each IP is symmetric with respect to 
the main direction of motion fixed by the vector w^x). Introducing the family of 
circular sectors 

S{x,R,a) = |y € : \y - x\ < R, f(x, y) ■ w{x) > cos || , 

where f(x, y) and 'w{x) are the unit vectors in the directions oi y ~ x and w{x), 
respectively, and setting 

i?c = max{i? > : ^ln{S{x, R, a^)) < p, C^i^Six, R, a^)) < s}, 

we can easily define B^i^x) :— S{x, Rc, ac). Notice that the constraint on the 
Lebesgue measure of each S basically amounts to fixing a metric upper bound 
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jjmax maximum radius allowed for the zone of cohesion Bc{x). This formal- 

izes the topological cohesion complemented with a metric cut-off observed in [3] 
(see also [TO). 

Conversely, we take the strength of the repulsive velocity v"^ proportional to the 
inverse of the distance among IPs within the zone of repulsion Br{x). Repulsion 
is mainly metric, as each IP simply tries to maintain a minimum distance between 
itself and other IPs. Hence Br{x) has a fixed size C'^{B,.{x)) , while its measure 
^n{Br{x)) may vary from point to point according to the crowding of the space. 
Assuming again for convenience that Br {x) is a circular sector of the ball with center 
in X and (fixed) radius Rr > 0, we simply have Br{x) :— S{x, Rr, otr), ar S (0, 27r] 
being the angular span of the zone of repulsion. 

Finally, the external velocity w depends essentially on the interactions of the IPs 
with the environment. Instead, it is independent of the distribution of the agents 
in the domain, because it represents the velocity that an isolated agent would set 
to reach its targets. Due to the strong differences in the structure of the typical 
environments where pedestrians and animals move, we refrain from giving here a 
general structure of w. We will detail it in the next subsections, with reference to 
specific applications. 

4.1. Macroscopic models. Macroscopic models are useful to study emergent self- 
organizing behaviors on large scales. This is the case of human crowds, in which 
self-organization entails the formation of patterns that are clearly visible only when 
the density of people is sufficiently high. In addition, these models are particularly 
handy to face the technical complexity of pedestrian flows in structured environ- 
ments scattered with obstacles, possibly also in connection with safety and opti- 
mization issues. 

Macroscopic models are obtained from the modeling framework outlined above 
via the continuum hypothesis, which amounts to saying that the measure /i„ is 
absolutely continuous with respect to the Lebesgue measure iHn < C"^)- There- 
fore, one can introduce a non-negative function pn G L^{il) such that dfin = Pn dx, 
and speak of density of IPs. The link with the measure /i„ is explicitly stated as 



E 

Starting from an initial measure ^ existence of the density for all times n 
is provided by the following result proved in [43 : 

Theorem 1. For all n > 0, let a constant Cn > exist such that 



If in addition po E L^{il) L°^{il) then pn E L^{il.) O L°"{il) as well, with 
1 1 Pn 1 1 oo — (niLi C'fe) IIpoIIoo for all n > 0. 



Under the continuum hypothesis, the intelligent velocity Q specializes as 





If Mo ^ ^'^ is nonnegative, then there is a unique sequence {pn}n>i C L^{il), 
Pn > a.e. in H., such that /x„ <C with dp,n — PndC^ and \\pn\\i = \\pq\\i for all 
n>0. 



(6) 



l^n[Pn]{x) 



I 



Ec{y - x) pn{y) dy + 



I 



Fr 



\y-x\ 



y~x 



2Pn{y) dy- 



B,(x) 
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Figure 3. Left: the grid {E'-'jjCi in the domain ft, with the 
generic cell highlighted. Right: the coefficients of the numerical 
scheme ([t]) are determined as the (hyper)volumes of the overlap- 
ping parts of adjacent grid cells under the push- forward operated 
by the approximate motion mapping 7^. 



Concerning the external velocity, to be definite we consider the application to 
pedestrians. In this case, it is convenient to model w as a normalized potential 
flow w = Vu/|Vm|, which does not depend on the distribution pn of the people but 
only on the geometry of the domain (cf. Sect. |3|, including the possible presence 
of obstacles. The function u : ^ M is a scalar potential satisfying Laplace's 
equation Au = 0. Boundary conditions for this equation may be set to identify 
attractive (resp., repulsive) areas, such as doors or displays (resp., obstacle edges 
or perimeter walls). For instance, targets along the boundary d^l of the domain 
may be characterized by the maximum potential, say u = 1, whereas the remaining 
portion of the perimeter walls by the minimum potential, say u — 0. On the 
internal boundaries, namely obstacle walls, the Neumann condition §^ = may be 
prescribed instead, which corresponds to zero normal component of the velocity w. 
Finally, the external velocity is generated by solving for u the resulting stand-alone 
elliptic problem 

Equations are numerically solved on a finite volume-like partition of the domain 
n, using an ad hoc computational scheme deduced from the one suggested in gi] 
with the inclusion of the new modeling features introduced here. Basically, a grid 
{Ej}j'Li consisting of elements of characteristic size /i > is introduced in 
(see Fig. [s] left), and the density /9„ and the velocity Vn are discretized as piecewise 
constant functions: 



Vx e E'l 



being a point of the cell ijj*, for instance its center. This induces naturally a 
discretization of the motion mappings 7„: 

ln{x)^tn{x)^X + vl[~p':^]{x)M, 
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where it should be noticed that the 7„'s act as piecewise translations on the grid 
{Ej}j^^. After defining the new measures dfl'^ :— p'^dC^, the scheme is ob- 
tained by imposing the push-forward fi^+i — In^P'n the grid cells E'j only, 
i.e., = /2^((7:^)-H4)) for all j = 1, . . . , M^ which gives (cf. Fig. | 

right) : 

If the 7„'s are sufficiently smooth, and if the spatial resolution h and the time step 
At are linked by a CFL-like condition, this scheme turns out to be nicely behaved 
in terms of stability and localization error: 

Theorem 2. Let 7„ be a diffeomorphism and let h, At > satisfy 

At max \vn[p'^]{x'^)\2 < h, 

] = l,...,Mh 

where | • I2 is the Euclidean norm in . Then: 

(i) There exists a constant C > independent of h such that 

A/'' 

J2 lMn+l(i?,") - ~^'n + liE^)\ < C {\\pn ^ P% + h) • 
J = l 

(ii) For each n > 0, there exists a constant Cn > independent of h such that 
max \iXn+i{E'^) - ~^':,+i{E';)\ < max 1^0(^5,') - 4(^)1 + 0,^. 

j — l,...,AI^ j — 

The proof of this result can be recovered again in [13]. 

The estimates of Theorem [2] rely essentially on the metric, because they 
assume the existence of the densities {pn}n>i for the measures {nn}n>i- However, 
it can be questioned that such a metric is not the optimal one to evaluate the 
distance between measures, even when densities are available. This is particularly 
true in the application to pedestrian flows, where self-organization phenomena may 
give rise to measures concentrated in thin areas of the domain (though not singular 
in view of Theorem [ij . Small errors in the localization of these measures would be 
roughly estimated by the distance of the corresponding densities, even when the 
corresponding distributions of the crowd are intuitively close. 

To be definite, let us consider two measures ni, p,2 with densities pi, p2, respec- 
tively, as illustrated in Fig. [4] This kind of distributions may be seen, for in- 



stance, along the vertical cross section of crossing pedestrian flows (cf. Sect. 5.1.1 1 
when oppositely walking lanes have fully emerged. Assume that ||/0i||oo = 0{l/e), 
£'*(supppi) = 0(e), i = 1, 2, for some arbitrarily small e > 0, whereas the average 
pointwise distance between pi and p2 is 0{S) for S > 0. Thus it is reasonable to 
expect a localization error of p2 over pi of the order of S, but if i5 is larger than e 
and we use the norm we inevitably get \\p2 — pi\\i = 0(1), because most mass 
is concentrated in the non-overlapping parts of the supports of pi, p2. 

A more correct way to measure the distance between pi, p2, which also better 
matches our intuition on what this distance should be, is offered by the Wasserstein 
metric, that we can briefly introduce as follows. Let (fi, d) be a metric space for 
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0(e) 



0(l/e) 




0{S) 

Figure 4. Two measures dfii = pidC^, i = 1, 2, whose dis- 
tance is 0(1) whereas Wasserstein distance is 0{5). The latter 
gives then a better estimate of the error produced by /i2 in the 
locahzation of /ii. 



which every probabihty measure on 51 is a Radon measure. Then the Wasserstein 
distance between two probabihty measures ni, ^2 is defined as 

^2) := inf j d{x, T{x))d^i{x), 
n 

the inf being taken among all transport maps from fii to fi2 '■ 

M2-T#Aii, I.e., ^i2iE) ^ fiiiT-\E)), y E e BiE). 

W{iii, 112) is the best (i.e., the lowest) transportation cost to move the measure 
111 onto /i2 (and vice- versa). In our previous example, we would actually get 
Willi, 112) = 0{5). 

The discrete-time model ([3| can be interpreted as an explicit Euler discretiza- 
tion in time of a gradient flow on the Wasserstein space, at least for sufficiently 
smooth motion mappings, cf. [U [SS] [SD] . Therefore the Wasserstein metric may be 
profitably used to further improve the theory with a more accurate error analysis. 
In particular, in future work we will try to obtain error estimates in space with 
respect to the Wasserstein distance, estimates in time being instead deducible from 

HI- 



4.2. Microscopic models. Microscopic models are useful to study self-organization 
phenomena at small scale. They are common in biological literature, where a wealth 
of models were studied in order to understand grouping behavior in fish schools, 
bird flocks, mammals herds and bacteria aggregations. Microscopic models can 
be obtained from our time-evolving measures framework, removing the continuum 
hypothesis. Denoting by a;" € the position of the j-th IP at time n, the measure 
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Figure 5. At the microscopic scale, the body size £ of the IPs 
matters. Therefore, repulsion from the j-th agent is always active 
in a complete neighborhood of radius £, regardless of the radius Rr 
and the angle a^. This, however, does not affect the zone of cohe- 
sion. 

/i„ is now chosen to be the counting measure, i.e., 

N 

(8) A^n = E ^^7 ' " ^ 

j=i 

where A'^ is the total number of agents and is the Dirac measure centered in x. 
Existence issues are now by far much easier than in the macroscopic case, because 
the measure fin+i can be constructed explicitly. Indeed, after a push-forward by 
7„ it is easily computed that 

N 

hence at the next time step the new positions of the IPs are x^^^ = ln{x'j). Finally, 
the resulting model turns out to be a classical agent-based model of the form 

which corresponds to the discrete-time version of the dynamical system Xj(t) — 
v(t, Xj(t)) obtained by an explicit Euler scheme. 

Inserting the measure ([s]) in Eq. Q yields the following expression of the intel- 
ligent velocity: 

(9) E ^c(x,"-x7)-f E ^^dt—^- 

Besides the formal derivation of Eq. ([9]) from Eq. (Q, we notice that at the 
microscopic scale the body size of the IPs matters. Indeed, the latter are not point 
agents, thus when two of them are closer than a characteristic distance > they 
touch each other (cf. also [S3 HHl HI] about microscopic pedestrians). Therefore, 
we admit that within a small neighborhood of radius (. repulsion from a certain 
agent is active against all other agents, regardless of the angular span a^, and 
that the distance between two IPs is never smaller than t, cf. Fig. [5j This also 
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helps avoid singularities, i.e., non-integrability with respect to the Dirac measure, 
in the second sum of Eq. ^ . 

Finally, we assume that the external velocity of the agents is constant (for in- 
stance, the vector w = (1, 0) for a rightward motion), as the environment is clear 
and we are mainly interested in the dynamics of the interactions among the IPs. 
This way, we can formally drop w from Eq. (|4]) by a simple change of frame of 
reference. 

5. Numerical results 

In this section we address some relevant case studies, which highlight the ability 
of our model to reproduce self-organizing patterns at both the macroscopic and 
the microscopic scale. In particular, we focus on pedestrians for macroscopic self- 
organization, studying the emergence of oppositely walking lanes in crossing flows, 
the spontaneous arrangement of people in a group in motion, and finally the dy- 
namics in a crowded environment scattered with obstacles. Conversely, we study 
microscopic self-organization with specific reference to animals, for which a closer 
look is necessary in order to catch crystal-like structures and line formations. 

5.1. Macroscopic self-organization in pedestrians. As recalled in Sect, [l] 
pedestrians experience, in normal situations, a very mild cohesion among each 
other compared to repulsion. Hence, in most numerical tests below we set Fc — 
and consider the following overall velocity: 

(10) ^"[^"K-) = ^+ J F.-l—^p.(y)dy, 

where < is constant. Due to the ahead-behind asymmetry of pedestrians, their 
visual field covers a frontal area only. Therefore, we define the zone of repulsion 
Br{x) to be the half-ball of radius Rr in the direction of the external velocity. 
Unless otherwise stated, the values of the relevant parameters are, for all tests, 

ar = IT, Rr = 0.1, Fr = —1. 

5.1.1. Lane formation in crossing flows. We consider two groups of pedestrians, of 
density pi^^^ e L^{fl), k = 1,2, respectively, walking in opposite directions. This 

(k) (k) 9 

amounts to defining two measures = pn dC , each of which evolves in time 

according to the push- forward ([s]). However, the two evolutions are coupled, since 
we assume that, within each group, repulsion is oriented against the individuals of 
the opposite group. Specifically, this means that the intelligent velocities are such 
that: 

-i^Hpi^^](^)= / Fr-^,p'^\y)dy, 
si) 

-i^npi^^](^)= / Fr-^,p^^\y)dy, 
J \y — x\ 

SO that pedestrians try to steer clear of oppositely walking people and to gain instead 
room in their direction. The external velocities are the vectors ui'^-* = (1, 0) and 
w^^^ = (—1, 0), thus groups walk rightward and leftward, respectively. 

Figure [6] shows that the model is able to account for lane formation when the 
two groups meet at the center of the domain and start to interact. In particular, 
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Figure 6. Lane formation in crossing flows, with inhoniogeneous 
uninterrupted (flrst row) and time-periodic (second row) injection 
of people from the boundaries of the domain. Pedestrians walking 
rightward are in blue, those walking leftward in red. Negative 
values of the density of the flrst ones are for graphical purposes 
only. 




Figure 7. Spontaneous arrangement of a crowd in motion. Start- 
ing from a compact cluster, the group expands and the density 
decreases. Leaders tend however to maintain the initial conflgura- 
tion. 



alternate lanes fully emerge for both uninterrupted and time-periodic flow of people 
from the boundaries of the domain, and turn out to be a quite stable equilibrium 
conflguration of the system. Such a conflguration is however reached in different 
times, being in particular more delayed in the second case, as a consequence of 
different intermediate dynamics undertaken by the system. Notice also the sponta- 
neous breaking of symmetry occurring between the two groups, which are instead 
specular at the beginning. 

5.1.2. Spontaneous arrangement of a crowd in motion. Next we investigate the self- 
organization spontaneously emerging within a group of pedestrians in motion in a 
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Figure 8. Effect of cohesion in walking pedestrians, (a) Initial 
condition, (b) Topological cohesion makes the three clusters merge, 
(c) Metric cohesion with large radius R"^^^ gives a result qualita- 
tively similar to the topological one. (d) Metric cohesion with small 
radius R^^^ produces a merging only of sufficiently close clusters. 



clear environment. The total velocity u„ is now as in Eq. (10 1, with an external 
velocity simply given by w = (1,0) (the group is walking rightward). The group 
is initially compact and has a homogeneous density (Fig. |7]i). As soon as people 
start to interact, the model predicts an expansion of the crowd in consequence of the 
repulsion. At the same time, the density decreases and becomes inhomogeneous due 
to the anisotropy of the visual field (Fig. ^jp). In particular, given the orientation of 
the external velocity, top-bottom symmetry is preserved, because _B,. is symmetric 
with respect to w. However, front-rear symmetry is lost, and most people remain 
initially concentrated in the rear part of the group, where the influence of the 
mass ahead is stronger. By consequence, in this zone the velocity is lower, hence at 
successive times the group elongates in the horizontal direction until the distribution 
of people becomes again substantially homogeneous (Fig. jTc). Only the motion of 
the leaders seems to be basically unperturbed (Figs. |7|d, [Tp), coherently with the 
fact that they simply follow the external velocity because nobody is in front of 
them. 

5.1.3. Effect of cohesion. In this test we investigate the effect of cohesion in the 
macroscopic framework of pedestrians. The intelligent velocity is now as in Eq. 
([6|, with Fc significantly greater than Fr {\Fc/Fr\ = O(IO^)) and a sensing domain 
for cohesion spanning the whole space around the agents (Oc = 27r). 

The test starts with three clusters of pedestrians at the same homogeneous den- 
sity, located a certain distance away from one another along the left side of the 
domain (Fig. [8^), which walk rightward {w — (1, 0)). If W^^^ is large enough, 
pedestrians are allowed to adjust the amplitude of their zone of cohesion Be so as 
to interact with a predefined amount of people, which in this simulation is set at 
I of the total mass initially present in the domain {i.e., p — |/io(f^) in Eq. (|5|). 
Then the three clusters tend to merge in a unique group, as shown in Fig. [8|d (this 
result can be compared with Fig. 4c in which shows the outcome of a similar 
experiment performed by a microscopic model). If instead p = +oo, the zone of 
cohesion is fixed to its maximum size determined by R^^^ (metric cohesion). In 
particular, for a large radius i?Jf (Fig. [s]:) the result is qualitatively similar to 
that obtained with topological cohesion, while for a small radius R™^^ (Fig- [sji) 
only the two clusters initially sufficiently close merge, the third one being instead 
unaffected by the presence of other agents in the domain (see Fig. 4b in [5]). We 
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(e) (f) (g) 

Figure 9. Pedestrian dynamics in presence of obstacles. Upper 
row: a crowd wants to reach the right edge of the domain, going 
through a botteneck. An obstruction forms as pedestrians try to 
access the passage, until all people flow to the opposite side. Lower 
row: speed map. The speed of the crowd is inhomogeneous, with 
a sharp transition from low to high values across the bottleneck in 
correspondence of the opposite transition in the values of the den- 
sity. 



refer the reader to the next Sect. |5.3| for comments on the importance of topological 
cohesion in spite of some qualitatively similar metric outcomes. 

5.1.4. Dynamics in presence of obstacles. Finally we study the motion of a crowd 
in a structured environment, in which some obstacles give rise to bottlenecks and 
direct pedestrians along preferential paths (the external velocity field is no longer 
homogeneous in space). In particular, we consider the case of a group of pedestrians 
wanting to go through a narrow passage, obtained by placing two obstacles in front 
of each other as in Fig. [9^. The external velocity w is obtained by solving Laplace's 
equation for the potential u, along with sliding boundary conditions at the obstacle 
edges (Neumann conditions). A Dirichlet boundary condition is instead imposed 
on the right edge of the domain, in order to set the potential at its maximum and 
to identify pedestrians' target. In this simulation cohesion is not active {Fc = 0), 
while repulsion is felt more strongly than in case of motion in clear environments 

{Fr = -10). 

Starting from an inhomogeneous density of people, confined in the left area of 
the domain ft (Fig. [9^), the self-organization of the crowd predicted by the model 
is as follows. When approaching the bottleneck (Fig. |9|d), people initially pass 
through at the maximum speed (Fig. |9^), however not all pedestrians can access 
the bottleneck at the same time and an obstruction forms (Fig. |9|:). Speed before 
the bottleneck is low, some individuals in the middle of the group are even forced to 
stop, whereas behind the bottleneck it attains again its maximum (Fig. [9j^). After 
a certain time, the whole group flows through the bottleneck and the obstruction 
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Figure 10. Two-dimensional globular clusters, metric interaction. 




(a) (b) (c) 

Figure 11. Two-dimensional globular clusters, topological inter- 
action, (a) Typical outcome with crystal structure, (b) angle dis- 
tribution among the agents for 1 run, and (c) angle distribution 
among the agents for 100 runs. 



is depleted (Figs. |9|i, [9^). This simulation compares qualitatively well with that 
proposed in [19J by means of a microscopic model. 

5.2. Microscopic self-organization in animals. For animal groups, we consider 
the complete structure of the intelligent velocity, in which both cohesion and re- 
pulsion are active. We investigate the importance of the topological correction and 
the effects of the anisotropic interaction, varying p and the angles and ac- We 
choose Rr as few times the body size £ of the agents, while we allow for a very large 
maximum radius i?™'^'^ for cohesion. 

5.2.1. Two-dimensional globular cluster. For the first test we set — — 2tt and 
p = N. The choice p = N, together with a large R^^^, implies that cohesion is 
basically metric and all-to-all. The parameters Fc and Fr are of the same order of 
magnitude. The system reaches a stable equilibrium in few iterations, forming a 
ball-like group with an irregular internal structure. In Fig. [lO]we show the typical 
outcome. If we introduce the topological correction, choosing p — 7, the system 
reaches again a stable equilibrium in few iterations, but this time a group forms, 
in which all IPs are at the same distance from each other. Every internal IP is 
surrounded by six group mates forming an hexagon, in a crystal-like structure (Fig. 
[ll}i) (cf also [161136]). 

By computing the distribution of the angles between each IP and its neighbors, 
we investigate the orientation of the hexagons. In Figs. [TT]d,[TT}: we show the angle 
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(a) (b) 

Figure 12. Macroscopic and microscopic models can reproduce 
the same patterns, confirming that they come from the same mod- 
ehng framework. 













(a) 



(b) 



(c) 



Figure 13. Line formations: (a) purely metric interaction, (b) 
topological interaction, and (c) angle distribution among the agents 
for 100 runs. 



distribution for 1 run and 100 runs, whence we can deduce that the orientation of 
the hexagons is mainly random. 

By switching cohesion off (Fc = 0) and setting up a mild repulsion among the 
agents (Fj. = —0.05) with a frontal visual field only (a^ = tt), we can also mimic with 
the microscopic model the spontaneous arrangement of a group of IPs described in 
Sect. 5.1.2| by the macroscopic model. Starting from a regular square configuration 
(Fig. 12[i), the IPs move rightward and interact only with the IPs in front of them. 
In Fig. |12|3 we show the final stable configuration, directly comparable with that 
in Fig. |7j:. 

5.2.2. Line formation. Here we set ar — j, ac — tt, and p = N. Cohesion is 
now greater than repulsion (Fc > \Fr\). The system does not reach an equilibrium 



(see Fig. 13 i), nevertheless the resulting pattern is strongly different from that of 



the previous test. The new outcome is mainly due to the modification in the in 
the relative magnitude of Fc and F^ and to different angles ac, ar, which imply a 
different anisotropy in the sensing zones. 

When introducing the topological correction with p — 7, the system reaches an 
equilibrium after few hundreds of iterations, forming a line oriented in the direc- 
tion of the motion. Neglecting small contributions of the repulsion component we 
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obtain a rather stable line (but for some little border effect in the head due to 



that group leaders cannot interact with p group mates ahead), see Fig. 13d. Lines 
are an example of pattern produced by self-organization of terrestrial animals like 
migrating penguins or elephants. 

A statistic on the distribution of the angles among the agents for 100 runs (Fig. 



13:) shows that a line configuration is always reached. 

In 13J the authors show the effect of a continuous variation of the pair (a^, etc) 
from (27r, 2tt) to (^, tt). As we have seen in the previous test, the first choice 
corresponds to a two-dimensional globular cluster. They find that wide angles 
induce a stretching of the group along the vertical direction, although in most runs 
the system does not reach an equilibrium. Conversely, small angles lead to a strong 
elongation of the group in the horizontal direction, and, in the limit case, to the 
formation of a line. 

5.3. Comments on the effect of topological correction. The previous tests 
clearly show that topological cohesion between IPs greatly changes the resulting 
pattern. This does not mean that it is impossible to obtain similar structures 
with a purely metric cohesion, by duly tuning the radius of cohesion (see e.g., the 



macroscopic test in Sect. 5.1.3). However, the topological correction is essential 
in order to deal with a large value of the maximum radius allowed. Indeed, as we 
have recalled in Sects. [2] [4] a metric upper bound R™^'^ to the radius of cohesion 
Rc exists, which translates the fact that IPs are in no case concerned with very 
far mates, and which should necessarily coincide with the fixed radius of cohesion 
in a purely metric approach. Now, R^^^ is in general rather large, because IPs 
are able to see quite far, and can be attracted even by far fellows if necessary. 
By consequence, once the group is formed, a purely metric cohesion with a large 
^max -^Q^j^jj imply attraction with an unreasonable number of other IPs, instead of 
feeling comfortable with the proper amount of IPs in the surroundings. Thus the 
topological correction is the only way to stay cohesive with a reasonable number of 
group mates while keeping a large maximum radius i?™^'^. 

As a further confirmation of this, the test on the effect of cohesion in the macro- 



scopic model (cf. Sect. 5.1.3 1 shows that a small value of R"^^^ in the purely metric 
approach distorts the cohesion itself: aggregation is only partial, since little far 
group mates may not be seen. 

6. Conclusions and research perspectives 

In this paper we have introduced a modeling framework for self-organizing intel- 
ligent particles, which takes into account both macroscopic and microscopic points 
of view, and is suitable to include topological and anisotropic interactions in an 
easy way and in any dimension. 

The results we have obtained from our numerical simulations suggest that these 
two features alone let self-organization emerge spontaneously, without forcing pat- 
tern formation via ad hoc agent-specific behavioral rules. In this respect, the most 
important parameters of the model are the ratio \Fc/Fr\, i.e., the relative strength 
of cohesive vs. repulsive terms, and the angles ac, Q^r, in other words the span of 
the sensing zones for cohesion and repulsion, respectively. 

The modeling technique by time-evolving measures that we have introduced is 
promising and deserves further investigation, because it enables one to address 
macroscopic and microscopic modeling by common mathematical structures and 
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tools. Specifically, we have used the macroscopic scale to study self-organization in 
pedestrian flows, as in this case clearly distinguishable patterns emerge only when 
the density of people is sufficiently large. In addition, the macroscopic approach 
may be profitably used in control and optimization problems connected to the im- 
provement of the flow and the safety of crowds. For instance, as recalled in Sect, [l] 
about Braess' paradox, this might imply optimization of the locations of some ob- 
stacles in crowded environments, like train stations or shopping malls. Conversely, 
we have studied self-organization in animals at the microscopic scale in order to 
catch the fine internal structure of the group and to highlight the appearance of 
regular structures (e.g., crystals) formed by few agents. The microscopic approach 
may be used to address problems in which the granularity plays an essential role. 
Furthermore, it may allow to introduce in the behavior of the agents stochastic 
effects which are not suited to an averaged macroscopic framework. 
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